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SUMMARY 


OBJECTIVE 

The  objective  is  to  illustrate  a  rational  approach  for  modeling  and  characterizing 
porous  electrodes.  Although  we  selected  the  metal  hydride  system  for  this  purpose, 
such  an  approach  is  applicable  to  other  systems  as  well  and  is  an  example  of  current 
trends  in  the  development  of  electrochemical  cells. 

RESULTS  AND  CONCLUSIONS 

With  an  increased  awareness  of  the  environment  and  how  human  activities  impact 
the  ecosystem,  current  trends  in  the  development  of  energy  storage  systems  are  to 
employ  less  hazardous  material.  Consequently,  recent  years  have  witnessed  increased 
interest  and  activity  in  the  development  of  the  “hydrogen  economy.”  There  hais  been 
particular  emphasis  in  using  hydrogen  as  the  fuel  in  energy  generating  devices  such  as 
hydrogen-oxygen  fuel  cells  or  as  rechargeable  batteries  in  which  the  negative  electrode 
contains  hydrogen  absorbed  in  an  alloy.  These  metal  hydride  systems  show  a  great 
deal  of  promise  due  to  their  high  theoretical  energj’  and  power  densities,  rechargeabil- 
ity,  and  potential  broad  range  applications.  Although  the  first  intermetallic  hydride 
was  developed  thirty  years  ago  and  despite  construction  of  the  first  prototype  metal 
hydrogen  electrode  twenty  years  ago,  metal  hydrogen  electrodes  still  experience  many 
shortcomings.  These  include  high  cost  of  alloys,  poor  hydrogen  storage  capabilities, 
difficult  activation,  pyrophoricity,  problems  of  impurities,  thermodynamic  instabili¬ 
ties,  and  corrosion  in  alkaline  media. 

Experience  has  shown  that  before  a  battery  system  is  ready  for  full-scale  appli¬ 
cations,  a  logical  progression  of  transitions  through  the  R&D  cycle  must  be  pursued. 
In  particular,  modeling  contributes  significantly  to  understanding  battery  operation. 
The  premise  of  modeling  is  that  the  description  of  electrode  operation  can  be  reduced 
to  a  set  of  precise  mathematical  statements  whose  solution  serves  as  a  guide  for  the 
selection  and  design  procedures  to  meet  performance  requirements.  A  blueprint  for  an 
effective  Ri:D  program  is  presented  which  (1)  emphasizes  a  rational  transition  from 
6.1  to  6.2  and  (2)  relies  on  the  synergistic  effects  arising  from  merging  theoretical 
concepts  with  the  requirements  of  practical  devices. 
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1.  INTRODUCTION 


With  an  increased  awareness  of  the  environment  and  how  human  activities  impact 
the  ecosystem,  current  trends  in  the  development  of  energy  storage  systems  are  to 
employ  less  ha7drdous  materials.  Consequently,  recent  years  have  witnessed  increased 
interest  and  activity  in  the  development  of  the  “hydrogen  economy.”  There  has 
been  particular  emphasis  in  using  hydrogen  as  the  fuel  in  energy  generating  devices 
such  as  hydrogen-oxygen  fuel  cells  or  as  rechargeable  batteries  in  which  the  negative 
electrode  contains  hydrogen  absorbed  in  an  alloy.  In  June  1992.  at  Argonne  National 
Laboratory,  DOE  sponsored  a  workshop  to  develop  a  research  strategy  for  a  new 
program  entitled  “Advanced  Battery  Technology  Research  and  Development”  (US 
Department  of  Energy,  1992).  This  workshop  focused  on  batteries  for  the  consumer 
market.  A  notable  feature  of  the  discussions  was  the  emphasis  placed  on  the  area 
of  rechargeable  batteries.  Among  the  many  systems  discussed,  rechargeable  batteries 
based  on  zinc,  lithium,  and  metal  hydride  negative  electrodes  and  nickel-oxide  positive 
electrodes  were  considered  leading  candidates  to  command  a  large  share  of  the  battery 
market. 

Experience  has  shown  that  before  a  battery  system  is  ready  for  full-scale  appli¬ 
cations,  a  logical  progression  of  transitions  through  the  R&D  cycle  must  be  pursued 
(Szpak  and  Smith,  1990).  In  organizing  this  report,  we  emphasize  the  role  of  model¬ 
ing  since  it  contributes  significantly  to  understanding  battery  operation.  The  premise 
of  modeling  is  that  the  description  of  electrode  operation  can  be  reduced  to  a  set  of 
precise  mathematical  statements  whose  solution  serves  as  a  guide  for  the  selection  and 
design  procedures  to  meet  performance  requirements.  The  objective  of  this  report  is 
to  provide  a  guide  for  the  development  and  optimization  of  the  performance  envelope 
of  a  battery  for  a  given  application  and  emphasizes  a  rational  transition  from  6.1  to 
6.2.  A  blueprint  for  an  effective  R  &  D  program  is  presented  which  relies  on  the 
synergistic  effects  arising  from  merging  theoretical  concepts  with  the  requirements  of 
practical  devices. 

2.  CURRENT  RESEARCH  IN  ELECTRODE 

DEVELOPMENT 

Nickel/metal  hydride,  NifMH,  sj^stems  show  a  great  deal  of  promise  due  to  their 
high  theoretical  energy  and  power  densities,  rechargeability.  and  potential  broad  range 
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applications.  Research  activities  in  the  construction  and  characterization  of  metal  hy¬ 
drogen  electrodes  were  discussed  at  the  International  Symposium  on  Metal-Hydrogen 
Systems  in  Uppsala,  Sweden  (June  8-12,  1992).  Although  the  first  intermetallic  hy¬ 
dride  was  developed  thirty  years  ago  and  despite  construction  of  the  first  prototype 
metal  hydrogen  electrode  twenty  years  ago,  metal  hydrogen  electrodes  still  experience 
many  shortcomings.  These  include  high  cost  of  alloys,  poor  hydrogen  storage  capa¬ 
bilities,  difficult  activation,  pyrophoricity,  problems  of  impurities,  thermodynamic  in¬ 
stabilities,  and  corrosion  in  alkaline  media.  In  reviewing  the  progress,  or  lack  thereof, 
in  the  field  of  metal  hydrides,  Suda  and  Sandrock  (1992)  concluded  that  new  and 
radical  changes  in  thinking  and  approach  were  required  to  overcome  the  difficulties 
encountered. 

The  renewed  interest  in  NifMH  systems  is  also  evident  from  the  topics  presented 
at  the  October  1993  meeting  of  the  Electrochemical  Society  held  in  New  Orleans, 
Louisiana.  Several  papers  dealt  with  the  role  of  the  interphase  in  H  sorption.  For 
example,  Ohnishi  et  al.  (1993)  reported  that  the  electrochemical  behavior  of  the  MH 
electrode  is  dependent  upon  the  surface  characteristics  rather  than  bulk  properties 
while  the  effect  of  surface  composition  on  sorption  of  hydrogen  from  both  gaseous  and 
liquid  phases  was  discussed  by  Florino  et  al.  (Florino,  Konstadinis,  Opila.  k  Fang, 
1993).  The  later  noted  that  MgO  affected  H  sorption  but  could  find  no  correlation 
between  the  film  composition  and  bulk  stability.  Similarly,  Ciureanu  et  al.  (1993), 
employing  both  electrochemical  (cyclic  voltammetry  and  impedance)  and  physical 
methods  (X-ray  diffraction,  Auger  spectroscopy),  concluded  that  the  cell  behavior  was 
primarily  governed  by  the  interphase.  Desilvestro,  Reymond,  Larcin,  and  Ruetschi 
(1993)  in  collaboration  with  Prof.  Schlapbach,  a  noted  researcher  of  MH  systems, 
examined  the  morphological  and  structural  changes  of  MH  electrodes  by  SEM,  X-ray 
diffraction,  and  ESC  A  techniques.  Modeling  of  the  discharge  cycle;  i.e..  modeling  of 
the  hydrogen  supply  to  the  reaction  surface,  has  been  discussed  by  Viitanen  (1992, 
1993)  and  Luo  et  al.  (Luo,  Johnson,  Feldberg,  k  Reilly,  1993).  In  the  latter,  discharge 
kinetics  were  taken  into  account. 

Modeling  of  metal  hydride  electrodes  requires  the  formulation  of  mathematical 
expressions  that  describe  the  transport  processes  within  and  between  both  the  elec¬ 
trolyte  and  the  hydrogen  absorbing  material.  Data  must  be  acquired  which  will  aid 
in  the  construction  and  verification  of  the  model.  Finally,  factors  that  affect  hydrogen 
sorption  must  be  examined. 
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2.1. 


ELECTRODE  MATERIAL  -  AB^  AND  AB^  ALLOYS 


There  are  two  classes  of  metal  hydride  alloys  employed  as  negative  electrodes, 
the  AB2  {ZrNi2  where  Ni  can  be  substituted  by  V,  Cr,  or  M?])  and  AB5  {La.\i5) 
alloys.  Both  classes  of  alloys  include  many  additional  metal  components  which  im¬ 
prove  performance,  life,  and  self-discharge.  With  the  passage  of  time,  what  at  first 
was  a  simple  binary  alloy  electrode  has  now  become  a  complex,  multicomponent  alloy 
system. 

Although  the  AB2  compounds,  such  as  ZrNi2,  show  reversible  hydrogen  absorp¬ 
tion  which  is  larger  than  that  of  LaNi^,  they  exhibit  poorer  charge  and  mass  transfer 
kinetics  (Zeuttel,  1992).  For  the  absorption-desorption  of  hydrogen  from  the  gas 
phase  for  both  the  AB2  and  ABs  alloys,  the  plateau  pressures  fall  between  10“^  and 
0.5  atm  at  20®  C  with  capacities  of  1.0- 1.3  HfMe.  This  translates  into  electrochem¬ 
ical  capacities  in  the  range  of  350  -  400  mA/g  (Canet  et  al.,  1992).  The  general 
composition  of  ABs-type  hydrogen  absorbing  alloys  used  for  battery  applications  is 
Mm{Ni  —  Co  —  Mn  —  Fe)x,  with  4.5  <  x  <  4.8,  where  Mm  denotes  the  “mis¬ 
chmetal”  (a  mixture  of  lanthanide  metals).  Depending  upon  the  battery  application, 
these  ABs-type  alloys  are  modified  by  the  addition  of  Mo,  B,  Ta,  W,  and  Zr.  For 
use  as  electrode  materials  in  the  construction  of  batteries,  these  alloys  must  exhibit 
corrosion  resistance  and  high  reactivity  for  the  charge  transfer  reaction  in  alkaline 
solutions  as  well  as  high  hydrogen  absorption  capacities  (Ise  et  al.,  1992).  The  AB5 
compounds,  in  particular  LaNi$  and  LaCo^,  are  not  stable  in  alkaline  electrolytes. 
Substitution  of  small  amounts  of  either  Ca  or  Ni  for  Co  has  been  found  to  improve 
the  stability  of  LaCo^.  This  stability  has  been  attributed  to  the  production  of  a  mul¬ 
tiphase  material.  The  problem  of  instability  in  alkaline  media  has  also  been  solved 
by  making  the  x-component  of  the  ABx  alloy  greater  than  5.  Besides  providing  long 
term  stability,  this  latter  modification  assures  higher  catalytic  activity  and  the  al¬ 
loys  are  easier  to  prepare  (Kanda  et  al.,  1992).  Today,  quality  controlled  electrode 
materials  are  available  from  Gesellschaft  fuer  Elektrometallurgie  m.b.H.,  Sparte  In¬ 
novation,  Hoefener  Strasse.  D-8500  Nuernberg  80  (Friedrich,  1992).  The  availability 
of  “made  to  order”  materials  greatly  accelerates  construction  of  electrodes  and  model 
verification. 

2.2.  ELECTRODE  PREPARATION 

Preparation  of  the  electrode  structure  is,  empirically,  one  of  the  most  important 
tasks  associated  with  the  development  of  practical  batteries.  Typically,  the  metal 
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hydride  electrode  for  a  sealed  Ni/MH  battery  consists  of  a  multicomponent  powder 
(50  MmNis  particle  size)  encapsulated  with  a  porous  0.5  pm  Ni  layer  deposited 
by  electroless  plating.  These  particles  are  placed  onto  a  ‘•foamed"  Ni  matrix.  I  sing 
electrodes  of  this  design,  AA  size  batteries  with  1.0-1. '2  AH  at  0.2  C  rate  were  built. 
This  is  equivalent  to  60  Wh/kg  or  170  Wh/1.  Charge  retention  was  found  to  be  70^ 
after  28  days  of  storage  at  room  temperature  (Yunshi  et  al..  1992).  A  somewhat 
different  electrode  preparation  was  described  by  Kuriyama  et  al.  (1992).  They  used 
three  methods  of  preparing  electrodes  of  an  A/mA^a.sCoo.TA/o.g  alloy.  The  particle 
size  of  the  alloy  ranged  between  0.12  and  0.15  mm.  In  preparation  A.  the  alloy 
particles  were  coated,  by  an  unspecified  process,  with  20  wt%  Cu  and  10  wt  % 
peril uororesin.  The  alloy  w'as  mixed  with  20  wt%  and  10  wt  %  perfluororesin  to 
yield  preparation  B.  In  preparation  C,  the  alloy  particles  were  simply  mixed  with 
10  wt  %  perfluororesin.  Each  mixture  was  compacted  into  a  pellet  and  provided  with  a 
current  collector.  Impedance  spectra,  from  0.5  mHz  to  10  kHz,  were  recorded  for  each 
electrode  at  various  discharge  depths  and  temperatures  under  open  circuit  conditions 
after  20  and  440  cycles  for  preparation  A,  after  20  and  410  cycles  for  B.  and  after 
20  and  130  cycles  for  C.  Electrochemical  impedance  spectroscopy  provides  a  means 
to  study  deterioration  mechanisms  of  porous  electrodes.  The  results  of  Kuriyama  et 
al.  (1992)  indicated  that  the  deterioration  of  C  resulted  from  poor  contact  between 
the  current  collector  and  the  pellet  as  well  as  between  the  alloy  particles  themselves. 
Furthermore,  the  Cu  layer  on  the  alloy  particles  used  in  preparation  A  was  found  to 
be  effective  in  preserving  good  electrical  contact  in  the  electrode. 

3.  PHYSICAL  PROPERTIES  OF  METAL 
HYDRIDE  SYSTEMS 

3.1.  KINETICS  OF  SORPTION/DESORPTION  OF  HYDROGEN 

Another  important  issue  in  developing  MH  electrodes  concerns  the  rate  at  which 
an  AB2  or  AB^  alloy  can  absorb  hydrogen  generated  either  electrochemically  or  sup¬ 
plied  from  the  gas  phase.  At  first  glance,  the  reported  data  on  the  subject  reveal  a 
great  deal  of  discrepancy  which  can  be  attributed  to  both  the  complexity  of  the  re¬ 
action  kinetics  and  the  transfer  path.  This  interpretation  is  based  upon  the  fact  that 
the  rate  constants  for  absorption  and  desorption  of  hydrogen  are  normally  reported 
under  the  condition  of  constant  temperature.  Since  absorption-desorption  processes 
are  either  endo-  or  exothermic,  this  condition  is  difficult  to  satisfy. 


4 


3.1.1.  Hydrogen  Sorption  from  Gas  Phase 


The  majority  of  the  studies  dealing  with  the  absorption  of  hydrogen  have  been 
conducted  from  the  gas  phase.  Dai*izer  et  al.  (Dantzer,  Pons,  Guillot,  Cai,  1992) 
and  Cimato  et  al.  (Cimato,  Zhang,  Goudy,  i:  Singleton.  1992)  reported  that  LoAjs 
absorbs  hydrogen  slower  than  Z,a.V24.7A/o.3-  This  observation  was  confirmed  by  Uchito 
and  Hoshimoto  (1992),  who  attributed  the  difference  to  surface  activity  promoted  by 
the  presence  of  /4/-oxides.  Stauffer,  Ezzehar,  and  Dreysse  (1992)  reported  on  an 
interesting  observation  of  hydrogen  clustering  near  the  metallic  surface.  Their  results 
indicated  that  surface  activity  is  very  specific  as  shown  by  the  differences  between 
Ni  and  Pd  surfaces.  Besides  the  existence  of  higher  hydrogen  concentrations  at 
the  Pd-surface,  there  are  two  energetic  states  of  adsorbed  hydrogen.  This  latter 
observation  has  support  from  recent  voltammetric  results  obtained  for  the  Pdf H2O 
system  (Szpak,  Mosier-Boss,  Scharber,  &  Smith,  1992;  Szpak,  Mosier-Boss,  Gabriel, 
&  Smith,  1994). 

3.1.2.  Hydrogen  Sorption  in  Electrochemical  Systems 

Due  to  the  non-autonomous  character  of  the  electrode/electrolyte  interphase,  a 
substantially  greater  number  of  factors  affect  the  sorption  process  than  in  the  gas 
phase.  An  interphase  region  is  formed  whenever  an  electrode  is  in  contact  with  an 
electrolyte.  In  the  simplest  case,  the  interphase  region  takes  the  form  of  the  electrical 
double  layer.  In  the  more  complex  cases  and,  in  particular,  during  the  charge  transfer 
reaction,  it  consists  of  layers,  each  associated  with  a  participating  elementary  process. 
In  this  representation,  the  interphase  region  is  an  open  system  in  which  a  number 
of  consecutive  processes  takes  place,  of  which  the  slowest  one  determines  the  rate. 
These  processes  include  transport  of  the  reactants  from  the  bulk  to  the  electrode 
surface  by  diffusion,  adsorption  on  the  electrode  surface,  charge  transfer,  desorption 
of  the  reaction  products,  followed  by  transport  of  the  reaction  products  away  from 
the  electrode  surface.  In  a  discharging  battery,  these  same  processes  occur;  however, 
in  a  battery  the  electrons  ultimately  flow  into  an  external  circuit  where  the  electrical 
work  is  delivered. 

In  the  negative  electrode,  the  relevant  processes  during  discharge  of  the  Ni  —  MU 
battery  occur  in  a  multiphase  environment  -  gas.  liquid,  and  solid.  Due  to  the 
hydrogen  absorption  capabilities  of  the  alloy  making  up  the  negative  electrode,  the 
electrode  itself  is  often  a  multiphase  system.  Transport  across  such  an  interphase  is 
a  dynamic  set  of  consecutive  processes  as  illustrated  in  Figure  1.  The  central  rol?  of 
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Figure  1.  Schematic  representation  of  an  interphase  for  a  hydrogen  absorbing  metal 
(fluxes  included),  (a)  Adsorption  plane;  (t)  charge  transfer  plane;  (1)  lattice. 


the  interphase  in  transport  of  electrochemically  generated  hydrogen  into  the  electrode 
interior  has  been  discussed  recently  (Szpak,  Mosier-Boss,  Scharber,  &  Smith,  1992; 
Szpak,  Mosier-Boss,  Gabriel,  &  Smith,  1994)  and  the  non-autonomous  character  of 
the  interphase  was  stressed.  It  was  shown  that  the  interphase  is  an  active  element  and 
that  its  properties  are  determined  by  those  of  the  contacting  phases,  i.e.,  the  electrode 
interior  as  well  as  the  electrolyte.  The  role  of  the  interphase  may  vary  from  that  of  an 
inhibitor  to  an  accelerator  of  both  charge  transfer  and  molecular  transport.  (Ohnishi 
et  al.,  1993;  Ciureanu  et  al.,  1993;  Dantzer,  Pons,  Guillot,  &  Cai,  1992;  Cimato, 
Zhang,  Goudy,  &  Singleton,  1992).  This  concept,  although  not  explicitly  stated, 
was  implied  in  the  discussion  of  hydrogen  absorption  from  the  gas  pheise  (Stauffer, 
Ezzehar,  &  Dreysse,  1992)  where  it  was  shown  that  small  metal  alloy  clusters  absorb 
more  hydrogen  via  a  mechanism  that  becomes  inoperative  as  the  cluster  size  increases. 
The  interphase  may  be  modified  in  the  course  of  battery  operation,  especially  as  a 
result  of  cycling.  For  example,  H^tsegawa  et  al.  (1992)  observed  the  development  of 
an  interphase  consisting  of  a  transition  metal  layer  which,  upon  cycling,  improved 
cell  performance. 

3.2.  THERMODYNAMICS  OF  SORPTION 

Recently,  thermodynamic  and  kinetic  factors  in  cathodic  //-sorption  (absorption 
and  adsorption)  into  metals  a. id  its  relation  to  //-adsorption  and  surface  poison¬ 
ing  effects  were  discussed  by  Jerkiewicz  (1992)  and  Jerkiewicz  and  Conway  (1992). 
Evidently,  the  operating  driving  force  for  sorption,  which  is  the  chemical  potential 
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gradient  across  the  interphase  region,  must  be  clearly  defined. 

3.2.1.  Chemical  Potential  of  Absorbed  Hydrogen 

According  to  Wagner  (1944),  the  chemical  potential  of  dissolved  hydrogen  is  given 
by 

Ai  = +  [/?T /n  +  Ap  (1) 

where  n  =  H/Pd  is  the  atomic  ratio  and  Aa<  represents  the  deviation  from  the  ideal 
distribution  of  H  among  the  octahedral  sites.  In  Wagner’s  formulation,  the  excess 
potential,  A/i,  is  due  to  the  contribution  of  hydrogen  electrons.  Later,  Brodowsky 
(1972)  added  the  protonic  contribution  so  that 

A/i  =  AAi+  +  A/i_  (2) 

resulting  in  the  following  picture  to  emerge;  //-atoms  within  the  octahedral  positions 
are  electronically  screened  protons  with  the  screening  charge  occupying  the  broadened 
states  in  the  4d-band  and  the  delocalized  states  in  the  5s-band  of  the  host  lattice. 
While  A/X+  is  generated  by  vibrating  protons,  A//_  is  the  electronic  contribution, 
which  increases  sharply  with  an  increase  in  //-atom  concentration.  In  reality,  the 
situation  is  further  complicated  by  lattice  expansion.  Furthermore,  changing  the 
number  of  defects  by  generation  of  new  defects  arising  from  the  lattice  expansion, 
affects  the  chemical  potential  and  through  it,  the  transport  into  the  bulk,  the  so- 
called  Gorski  effect  (Voelkl,  1972). 

3.2.2.  Chemical  Potential  of  Hydrogen  in  the  Interphase 

The  chemical  potenti2il  of  molecules  residing  in  the  interphase  region  must  reflect 
its  non-autonomous  character.  Jerkiewicz  and  Conway  (1992)  treated  the  sites  of 
adsorbed  and  absorbed  H  in  terms  of  site  fraction  statistics  in  order  to  write  down 
their  chemical  potentials.  The  efficiency  of  absorption,  in  relation  to  the  hydrogen 
evolution  reaction,  can  be  ex2unined  in  terms  of  the  reaction  mechanism.  The  effect  of 
catalytic  poisons  promoting  the  absorption  is  discussed  in  terms  of  competition  for  the 
sites.  The  discussion  focused  on  the  complex  nature  of  the  interphase,  however,  it  did 
not  model  the  electrode  response  to  the  charging  (discharging)  currents.  Within  an 
electric  field  (i.e.,  the  electric  double  layer)  there  is  an  increaise  in  energy  in  the  system 
which  has  not  been  included.  Since  the  chemical  potential  of  a  molecule  is  affected  by 
the  electrode  polarization,  the  picture  presented  b)'  Jerkiewicz  and  Conway  (1992)  is 
incomplete  and  additional  work  is  needed.  In  summary,  the  predictive  capabilities  of  a 
model  are  determined  by  several  factors,  among  them,  the  structure  of  the  interphase 
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and  a  full  description  of  the  operating  driving  force,  including  a  complete  expression 
for  the  gradient  of  the  chemical  potential,  both  in  the  bulk  electrode  and  across  the 
interphase. 

4.  TECHNOLOGY  ISSUES 

The  development  and  construction  of  metal  hydride  electrodes  for  battery  ap¬ 
plications  encompasses  four  areas  of  research;  (i)  the  rate  at  which  hydrogen  can 
be  recovered  and  brought  into  the  reaction  zone,  (ii)  the  rate  of  charge  transfer, 
(iii)  the  determination  of  the  amount  of  hydrogen  stored  within  the  electrode  struc¬ 
ture,  and  (iv)  reversibility.  The  first  and  second  areas  of  research  are  concerned  with 
the  elementary  processes  and  modeling  while  the  third  deals  with  the  hydrogen  sorp¬ 
tion  capacity.  The  fourth  area  of  research  deals  with  the  degradation  of  the  electrode 
upon  cycling. 

4.1.  THEORETICAL  ASPECTS 

The  cell /electrode  modeling  follows  the  rules  first  defined  by  Newman  and  Tobias 
(1962)  and  somewhat  later  extended  to  complex  systems  exhibiting  unique  features 
by,  among  others,  Tsaur  and  Pollard  (1984a,  1984b).  A  summary  of  the  govern¬ 
ing  equations  for  the  Tsaur-Pollard  formulation  for  the  Lif  S0Cl2  system  is  given  in 
Table  I.  With  appropriate  modifications,  these  equations  can  be  adapted  to  represent 
the  jV/if-system.  The  resultant  modifications  must  reflect  the  characteristic  features 
of  the  Af  if -electrode  as  displayed  by  experimental  simulation.  Furthermore,  the  need 
to  include  a  synergistic  interaction  between  the  assumed  events  as  shown  from  mod¬ 
eling  and  the  behavior  observed  in  experimental  analogs  cannot  be  overemphaisized. 

4.1.1.  Review  of  MH  Electrode  Modeling 

To  date,  only  certain  aspects  of  MH  electrodes  have  been  modeled.  In  gen¬ 
eral,  these  modeling  approaches  have  treated  the  MH  electrode  as  a  planar  sur¬ 
face  in  contact  with  an  electrolyte  separated  by  an  interfacial  region.  For  example, 
Luo  et  al.  (Luo,  Johnson,  Feldberg,  &  Reilly,  1993)  used  a  model  based  on  diffu¬ 
sion  out  of  a  cylinder  to  describe  the  kinetics  of  discharge  of  the  Pdo.gsA’^io.is  hydride 
electrode.  In  a  separate  exercise,  Szpak  et  al.  (Szpak,  Gabriel,  Smith,  &  Nowak, 
1991)  derived  equations  governing  transport  of  H  across  the  interphase  into  the  bulk 
electrode  for  the  PdlH20  system  and  have  obtained  numerical  solutions.  In  their 
model,  they  took  into  account  the  processes  illustrated  in  Figure  1.  This  model 
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Table  I.  Summary  of  governing  equations  for  Lij SOCI2  cells.  (Adapted  from  Tsaur  and  Pollard,  1984a  and  1984b). 
List  of  symbols  and  reaction  equations  for  Li/SOCh  electroreduction  are  found  in  sections  6.1  and  6.2,  respectively. 
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Equation  to  describe  dSj  2Dj  ,  l  i  \ 

thickness  of  defective  ~  ^LiCl  ~7  (Co)j  —  (Ca) 

LiCl  film  at  anode  ^ 


assumed  no  geometrical  restrictions  on  the  electrolyte  phase  and  agreed  well  with  the 
experimental  data  obtained  by  Riley  et  al.  (Riley,  Seader,  Pershing,  k  Walling.  1992) 
as  was  demonstrated  by  Szpak  et  al.  (Szpak,  Mosier-Boss,  Gabriel,  k  Smith,  1994). 

There  has  been  some  effort  to  place  more  realistic  constraints  upon  the  system 
when  modeling.  Recently  Viitanen  (1992,  1993)  has  modeled  the  transport  of  H  in 
a  MH  electrode  comprised  of  spherical  particles  of  uniform  size.  Mao  and  White 
(1991)  modeled  the  self-discharge  of  a  NiOOH  electrode  in  a  hydrogen  environment. 
This  model  took  into  account  diffusion  of  dissolved  hydrogen  in  an  electrolyte  film 
which  c  v'ers  a  flooded  electrode,  electrochemical  oxidation  of  hydrogen,  reduction  of 
NiOOH,  and  changes  in  surface  area  and  porosity  of  the  electrode  during  the  ‘ 
discharge  process.  Their  model  does  not  specifically  pertain  to  the  MH  electro 
however,  it  does,  peripherally,  describe  transport  of  an  electrolye  phaise  in  a  porous 
structure. 


4.1.2.  Limitations  of  Current  Models 

Clearly,  previous  modeling  of  the  M H  system  has  been  somewhat  incomplete  and 
limited  to  specific  elements  of  the  charging/discharging  cycle.  To  our  knowledge, 
no  general  model  has  been  developed  to  describe  the  transport  processes  within  and 
between  the  electrolyte  and  the  hydrogen  absorbing  material  during  the  course  of 
cell/battery  operation.  Ideally,  such  a  model  would  describe  every  aspect  of  MH 
charging/discharging.  This  is  essential  in  order  to  predict  the  operational  charac¬ 
teristics  of  the  system,  as  a  whole,  and  to  assess  the  influence  of  changes  in  design 
p£irameters  on  the  cell  performance.  Not  only  would  such  a  model  identify  system 
limitations,  but  it  would  also  help  guide  future  experimental  research. 

As  summarized  in  Table  I,  when  modeling  a  charging/discharging  cell,  one  must 
take  into  account  the  material  balances,  Faraday’s  law,  Ohm’s  law,  and  the  polariza¬ 
tion  equation.  For  the  M H-ele*  trode,  we  also  have  to  include  the  energy  balance,  the 
effect  of  the  interphase,  and  H  transport  between  and  within  the  particles  making 
up  the  porous  structure  of  the  MH  electrode.  These  latter  effects  will  be  taken  into 
account  by  the  form  of  the  auxiliary  equations.  The  development  of  a  general  model 
is  a  nontrivial  undertaking  since  charging/discharging  of  a  metal  hydride  is  a  very 
complicated  process,  as  shown  in  Figure  1,  and  there  are  essentially  three  phases  to 
consider  -  the  electrolyte,  the  particles  making  up  the  M H  electrode,  and  bubble  for¬ 
mation.  Therefore,  the  appropriate  equations  would  have  to  be  formulated  to  include 
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coupling  of  the  electrode/electrolyte  processes  within  this  multiphase  environment. 
From  our  viewpoint,  previous  modeling  of  the  MH  system  suffers  from  a  number  of 
deficiencies.  They  are: 

(i)  conservation  of  energy  from  which  the  temperature  distribution  is 
determined 

(ii)  the  interphase  which  is  an  active  element  in  the  transport  of  H  and  is 
non-autonomous  (Ohnishi  et  al.,  1993;  Ciureanu  et  al.,  1993;  Dantzer, 
Pons,  Guillot,  k.  Cai,  1992;  Cimato,  Zhang,  Goudy,  k  Singleton,  1992). 

A  general  model  will  have  to  include  the  effect  of  a  changing  interphaise 
on  H  transport.  The  earlier  model  of  the  Pd! H2O  system,  by  Szpak 
et  al.  (Szpak,  Gabriel,  Smith,  k  Nowak,  1991)  partially  treated  this 
effect  but  will  have  to  be  extended. 

(iii)  the  effects  of  the  electrode  matrix  and  grain  boundaries  upon  H  trans¬ 
port.  The  electrode  matrix  refers  to  the  porous  electrode  structure 
while  the  grain  boundaries  refer  to  the  metallurgical  ^lspects  within 
individual  particles.  Consequently,  we  must  deal  with  H  transport 
between  and  within  particles  making  up  the  porous  structure.  Trans¬ 
port  within  grain  boundaries  has  been  partially  treated  by  Zaeschmar 
(1983).  His  formalism  will  require  modification  to  be  applicable  for  the 
MH  electrode. 


For  the  construction  of  the  auxiliary  equations,  experimental  data  are  needed  re¬ 
garding  electrolyte  transport  properties  and  their  dependence  on  temperature  and 
composition,  properties  of  the  metal  hydride  as  a  function  of  loading  (e.g.,  conduc¬ 
tivities,  densities,  diffusion  coefficient  of  H,  etc.),  kinetic  parameters  (e.g.,  exchange 
current  densities,  active  specific  surface  areas,  symmetry  factors,  and  morphology 
parameters),  etc.  Although  there  are  experimentzd  dischau-ge  curves  for  working  cells 
in  the  literature  and  a  wealth  of  information  on  the  kinetics  of  charge  transfer  for 
many  half-cell  reactions,  pertinent  input  data  on  transport  processes  within  porous 
structures  are  lacking. 

Intuitively,  one  would  expect  interactive  cooperation  between  modeling  exercises 
aind  experiment  since  model  construction  and  verification  require  that  the  calculated 
results  match  actual  experimental  data.  Ideally  such  an  interaction  between  model¬ 
ing  and  experiment  would  be  synergistic;  experiments  are  designed  and  conducted  to 
provide  the  data  needed  to  aid  in  formulating  the  auxiliary  equations  in  the  model. 
Similarly,  the  model  may  predict  a  certain  behavior  which  will  have  to  be  verified  by 
experiment.  Rarely  does  this  synergistic  interaction  between  modeling  ^lnd  experi¬ 
ment  occur.  Clearly,  this  is  another  area  in  modeling  which  is  deficient.  When  the 
experimental  data  are  not  available,  it  is  a  common  practice  in  modeling  to  employ 
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the  general  Newman-Tobias  formalism  (Newman  L  Tobias,  1962)  and  simply  vary  the 
parameters  in  such  a  way  that  the  results  from  the  model,  e.g.,  a  calculated  discharge 
curve,  approximates  that  measured  for  a  working  cell.  Dependent  upon  the  number 
of  parameters,  it  is  quite  possible  that  a  large  number  of  solutions  exist  which  would 
approximate  the  observed  cell  behavior. 

4.2.  SPECIFIC  TASKS  FOR  MODELING 

Developing  a  generad  model  of  the  charging/discharging  MH  electrode  would  re¬ 
quire  (i)  modification  of  the  equations  in  Table  I  to  describe  transport  processes  within 
and  between  both  the  electrolyte  and  the  hydrogen  absorbing  material; 
(ii)  designing  and  building  experimental  analogs  to  secure  such  data  as  necessary 
for  the  construction  and  verification  of  the  model;  and  (iii)  examination  of  the  factors 
affecting  hydrogen  sorption. 

4.2.1.  Development  of  Formalism  for  the  MH  Electrode 

In  developing  a  general  model  to  describe  the  A/ if -electrode,  we  start  with  the 
concept  that  the  processes  occurring  within  a  porous  electrode  correspond  to  those 
on  a  planar  electrode  subject  to  the  geometrical  constraints  imposed  by  the  structure 
itself.  A  one  dimensional  model,  consisting  of  compacted  spheroid  particles  in  contact 
with  a  current  collector  on  one  side  and  continuous  electrolyte  phase  on  the  other,  is 
assumed.  Initially  spherical  particles  of  uniform  size  are  considered.  However,  later 
modifications  of  the  model  would  have  to  take  into  account  particles  of  non-uniform 
size.  The  porous  regions  of  the  cell  are  assumed  to  be  flooded. 

The  fluxes  involved  in  charging  the  metal  hydride  electrochemically  are  summa¬ 
rized  in  Figure  1.  When  charging  the  MH  electrode,  reduction  of  water  must  first 
occur,  via  the  Volmer  path,  Equation  (I),  or  the  Heyrovsky-Horiuti  path,  Equation 

(II), 


H2O  +  e- 
H2O  +  H{a)  +  e- 


H{a)  +  OH-  ;  (I) 

-+  H{a)  +  OH--,  k2,k^2,V  (H) 


with  gas  evolution  occurring  by  the  Tafel  path,  Equation  (III), 
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H{a)  +  H{a)  -*  H2{s)  ;  (Ill) 

Here  (a)  and  (3)  denote  species  present  in  an  adsorbed  state  and  in  solution,  re¬ 
spectively.  The  absorption  into  the  lattice  of  electrochemically  generated  hydrogen  is 
given  by  Equation  (IV) 

(IV) 

where  (/)  denotes  the  species  in  the  lattice. 

For  discharge,  the  reactions  tahing  place  at  the  surface  of  the  metal  particles  are 
the  reverse  of  reactions  (IV)  and  (I).  The  reverse  of  reaction  (IV)  is  the  dehydridiza- 
tion  reaction,  where  the  hydrogen  is  released  from  the  hydride  phase,  resulting  in 
the  formation  of  adsorbed  hydrogen  atoms.  The  electrochemical  reaction,  in  which 
current  is  generated  when  the  adsorbed  hydrogen  reacts  with  hydroxide  ions,  is  the 
reverse  of  reaction  (I).  Gas  evolution  occurs  via  the  Tafel  path,  Equation  (III). 

Returning  to  the  summary  of  equations  presented  in  Table  I,  the  material  balance 
equation  for  the  electrolyte  species  requires  no  further  modification  for  the  M H  sys¬ 
tem.  However,  the  solid  phase  material  balance  equation  will  have  to  be  modified  to 
take  into  account  the  change  in  volume  of  the  particles  during  charging/discharging 
and  grain  size  boundaries.  During  charging/discharging,  the  metal  hydride  particles 
expand/contract  resulting  in  porosity  changes.  The  upper  and  lower  limits  would  be 
the  total  volume  of  particles  fully  loaded  with  hydrogen  and  virgin  material  which 
has  not  undergone  a  charge/discharge  cycle,  respectively.  Once  a  hydrogen  absorbing 
material  has  undergone  a  charging  cycle,  the  lattice  does  not  collapse  back  to  its 
pre-loaded  condition.  The  degree  of  expansion/contraction  of  the  particles  making 
up  the  MH  electrode  is  a  function  of  the  degree  of  loading.  Modifying  the  solid- 
phase  material  balance  to  take  into  account  grain  boundaries  is  more  difficult.  Grain 
boundaries  represent  locations  of  a  high  density  of  lattice  imperfections.  The  major 
effect  grain  boundaries  have  on  charging/discharging  is  on  the  diffusion  of  hydrogen 
into  the  lattice.  The  interaction  of  an  interstitial  with  a  lattice  defect  is  energetically 
quite  different  from  its  interaction  with  the  undisturbed  lattice.  One  approach  is  to 
treat  grain  boundaries  as  internal  surfaces  on  the  spherical  MH  particles  where  an 
exchange  process  similar  to  adsorption-desorption  can  occur. 

Although  the  flux  relation  equation  for  the  electrolyte  species  in  Table  I  can  be 
used  for  the  MH  electrode  without  further  modification,  simplifications  may  result 
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from  the  introduction  of  effective  parameters,  e.g.,  Z?e//-  However,  a  corresponding 
flux  relation  equation  will  have  to  be  derived  for  diffusion  of  hydrogen  into  the  solid 
state  particles.  This  latter  equation  will  be  dependent  upon  the  surface  coverage  of 
adsorbed  hydrogen  atoms. 

The  expression  for  Faraday’s  law  in  Table  I  is  applicable  to  the  M H  electrode. 
However,  the  expression  for  Ohm’s  law  will  have  to  be  reformulated  to  comply  with 
the  MH  system.  The  polarization  equation  in  Table  I  will  also  have  to  be  extended 
to  include  changes  in  the  reaction  path  resulting  from  the  influence  of  an  active  inter¬ 
phase.  The  auxiliary  equation  in  Table  I  which  describes  the  thickness  of  the  defective 
LiCl  film  at  the  anode  is  not  applicable  to  the  MH  system  and  will  be  excluded. 
Conversely,  the  auxiliary  equation,  which  describes  the  active  specific  surface  area  in 
the  cathode,  will  be  modified  for  the  MH  electrode.  In  particular,  the  specific  surface 
area,  a,  is  a  function  of  c,  c,,  and  j.  This  latter  equation  will  also  have  to  include  the 
effect  of  gas  evolution  since  gas  bubbles  will  effectively  block  the  active  sites  on  the 
particles. 

An  additional  equation  is  required  to  take  into  account  the  energy  balance.  A 
comprehensive  treatment  of  thermal  effects  in  electrochemical  cells  is  given  by,  among 
others,  Newman  (1973).  As  a  starting  point,  the  temperature  distribution  is  obtained 
by  assuming  a  conservation  of  energy.  The  state  of  any  system  is  completely  char¬ 
acterized  by  the  internal  energy,  U;  the  volume,  V;  and  the  mole  number,  nj.  For 
systems  under  constant  pressure,  the  use  of  the  enthaJpy  function  is  more  convenient. 
The  transition  from  energy  to  enthalpy  representation  is  made  by  substituting  volume 
with  pressure  so  that  17^“^  = 

Identifying  the  porous  structure  as  an  open  system  and  applying  energy  conser¬ 
vation  in  the  enthalpy  representation  to  a  subsystem,  q,  the  change  in  the  enthalpy 
is  given  by 

where  a  and  ^  designate  subsystems  in  contact  with  each  other.  In  equation  (3),  dq 
is  the  heat  of  reaction  (e.g.,  absorption),  dw  is  the  work  done  by  the  system  (e.g., 
electrical  work),  and  is  ihe  heat  transferred  from  subsystem  q  to  subsystem 

The  temperature  distribution,  T{x,i),  is  calculated  from  the  change  in  the  enthalpy 
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(4) 


occurring  within  time,  dt,  viz., 

^  -  dJi^^  +  (f 

where  a  denotes  the  position  in  co-ordinate  x.  By  identifying  subsystem  q  with  po¬ 
sition,  we  identify  a  continuous  system,  as  shown  in  Figure  2.  A  complete  treatment, 
including  mass  and  momentum  balance,  leading  to  the  temperature  distribution  is 
given  in,  e.g.,  Haase  (1969). 


porous  electrode  structure 


Figure  2.  One-dimensional  model  of  a  porous  structure  showing  Q,  subsystems. 


The  effect  of  an  active,  non-autonomous  interphaise  will  be  taken  into  account  by 
an  auxiliary  equation.  Inclusion  of  diffusional  transport  of  hydrogen  into  the  lattice 
through  an  active  interphase  necessitates  the  separation  of  the  interph«ise  region  from 
the  diffusional  space.  In  their  model  of  charging  of  the  Pd  lattice,  Szpak  et  al.  (Szpak, 
Gabriel,  Smith,  &  Nowak,  1991)  treated  the  bulk  electrode  as  a  body  consisting  of  N 
layers;  the  layer  having  a  volume  K  and  the  area  of  the  dividing  surface  between 
it  and  the  (e  -f-  1)*^  is  5,,,+i.  In  general,  Vi  ^  VJ+j  and  ^  The 

diffusion  of  hydrogen  in  the  bulk  of  the  solid  body  was  considered  to  be  the  process 
of  jumping  from  one  layer  to  an  adjacent  layer  with  a  rate  that  is  proportional  to 
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the  interfacial  area  between  the  layers,  as  well  as  the  number  of  vacant  sites  in  the 
final  layer.  The  boundary  conditions  are  determined  by  the  physical  considerations  of 
events  occurring  at  the  1**  and  layers.  This  treatment  will  have  to  be  extended  to 
an  interphase  which  actively  participates  in  hydrogen  transport  into  the  lattice  and 
whose  structure  changes  throughout  the  charging/discharging  cycle.  As  the  loading 
increases/ decreases  during  charging/discharging,  the  number  and  thicknesses  of  the 
layers,  N,  making  up  the  interphase  vary  as  does  the  number  of  available  sites  within 
a  layer  and  the  jump  rate  constant  between  layers. 

4.3.  EXPERIMENTAL  ASPECTS:  MODEL  CONSTRUCTION  AND 
VERIFICATION 

The  exact  form  of  the  expressions  used  to  describe  the  MH  electrode  during 
charging/discharging  will  be  dependent  upon  the  measured  response  of  experimental 
analogs.  These  experimental  analogs  are  used  to  simulate  the  behavior  of  a  working 
cell,  thereby  providing  the  input  data  required  for  construction  and  verification  of  the 
model.  The  input  data  needed  for  model  construction/verification  are: 

(i)  electrolyte  transport  properties  and  their  dependence  on  temperature 
and  composition 

(ii)  solid  phase  conductivities  and  densities 

(iii)  metallurgical  aspects  of  the  MH  electrode  as  a  function  of  loading 

(iv)  the  open-circuit  cell  potential  and  its  dependence  on  temperature 

(v)  hydrogen  absorption  capacities  and  rates  of  ch2urging/discharging 

(vi)  kinetic  parameters:  e.g.,  exchange  current  densities,  active  specific  sur¬ 
face  areas,  symmetry  factors,  morphology  parameters,  and  effects  of 
the  interphase 


With  this  information  at  hand,  the  model  can  be  used  to  predict  the  time  and  position- 
dependent  behavior  of  the  system. 

4.3.1.  Simulations 

For  realistic  representation,  the  proposed  model  must  be  based  on  analyses  of  the 
actual  elementary  processes  occurring  within  the  confines  of  the  porous  structure.  Ex¬ 
perimental  data  on  local  and  bulk  transport  and  charge  transfer  kinetics  are  needed. 
Similarly,  data  are  needed  on  the  variation  in  the  driving  force  which  occurs  natu¬ 
rally  as  a  function  of  distance  as  well  as  morphological  changes  for  any  given  mode 
of  operation.  Such  data  are  difficult  to  obtain  in-situ  since  the  insertion  of  probes 
compromises  the  integrity  of  the  porous  structure.  Instead  we  propose  to  use  simpler 
experimental  arrangements  that  are  designed  to  simulate  the  behavior  of  practical 
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electrodes.  The  intent  of  simulation  is  to  represent  the  real  system  by  one  which  is 
accessible  to  direct  observation  so  that  the  effects  of  varying  structural,  electrochemi¬ 
cal,  and  transport  phenomena  can  be  followed  using  ultramicroelectrodes  and  optical 
techniques.  In  the  simulation  attempt,  we  select  a  configuration  that  exhibits  some 
or  all  of  the  features  common  to  the  operating  system.  An  example  of  the  evolution 
of  an  experimental  analog  to  represent  a  single-pore  is  illustrated  in  Figure  3. 


anode  cathode  anode  cathode 


ANODE  CATHODE 


Figure  3.  Schematic  derivation  of  single-pore  analog  cell  through  successive  degrees 
of  abstraction:  A,  real  electrode  consisting  of  interconnected  pores  and  bonded  parti¬ 
cles;  B,  idealized  abstraction  of  A,  retaining  tortuous  pores,  but  pores  are  no  longer 
interconnected;  C,  further  abstraction  where  pores  pass  straight  through  the  electrodes 
(tortuosity  factor  =  1);  D,  section  through  the  cell,  shown  as  an  abstraction  and  as 
the  analog  single-pore  cell. 


Since  we  are  interested  in  obtaining  data  on  transport  processes  within  a  porous 
structure,  we  would  use  a  segmented  electrode  similar  to  that  shown  in  Figure  4  as 
an  experimental  tool.  This  particular  experimental  design  wzis  used  to  study  the 
current  density,  dependence  of  the  AgfAgCl  system  (Szpak  &:  Katan,  1975; 

Szpak,  Nedoluha,  &:  Katan,  1975).  In  the  construction  of  the  experimental  analog,  a 
silver  working  electrode  consisting  of  several  segments,  as  shown  in  Figure  3,  is  vapor 
deposited  onto  a  glass  microscope  slide.  The  layer  of  silver  will  act  as  a  current  col- 
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lector  and  a  substrate  onto  which  a  hydrogen  sorbing  material  will  be  deposited.  The 
silver  has  the  added  advantage  of  acting  as  a  barrier  to  diffusion  of  sorbed  hydrogen 
within  the  absorbate.  Since  there  is  a  wea’th  of  data  available  for  electrochemical 
loading/unloading  of  Pdf  H,  this  system  is  often  employed  as  a  model.  In  addition, 
some  commercially  available  MH  materials  for  battery  applications  are  coated  with  a 
catalytic  surface  of  Pd.  Therefore,  initially,  the  hydrogen  sorbing  material  would  con¬ 
sist  of  Pd,  either  vapor  or  electrochemically  deposited  onto  the  Ag  substrate.  Later 
commercially  available  AB2  and  ABs  materials  would  be  employed  once  the  experi¬ 
mental  details  have  been  established.  These  materials  are  particulate  and  will  have 
to  be  mixed  with  Ag  to  provide  a  conductive  path,  sintered,  and  connected  to  the 
segmented  Ag  current  collector  on  the  glass  slide.  In  order  to  make  direct  measure¬ 
ments  of  the  space-time  distributions  of  transfer  current,  ji(i,  t),  the  thin  segments  of 
the  electrode  structure  must  bti  isolated  and  provided  with  a  separate  current  take¬ 
off  lead.  This  should  be  done  so  that  the  sectioning  will  not  appreciably  change  the 
structure  and  elementary  processes  that  comprise  the  over-all  reaction  path.  Also 
provisions  must  be  made  to  aissure  the  monitoring  of  current  in  each  segment  without 
adtering  the  original  j{x,i)  distribution. 


Figure  4.  Experimental  cell  assembly.  Lower  right:  details  of  holder  construction, 
electrical  connections,  and  positioning  of  segmented  electrode.  Upper  left:  construction 
of  segmented  electrode  showing  vapor  deposited  silver  on  a  glass  substrate.  The  shaded 
area  is  a  U-shaped  Teflon  tape,  60  fim  thick. 


The  experimental  cell  assembly,  illustrated  in  Figure  4,  lends  itself  to  additional 
measurements  besides  that  of  These  include  both  electrochemical  as  well  as 
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optical  measurements.  For  example,  the  rate  of  H  absorption/desorption  can  be 
studied  under  both  potentiostatic  and  galvauostatic  conditions.  The  working  elec¬ 
trode  can  be  pulsed  either  potentiostatically  or  galvanostatically  while  the  behavior 
of  the  segments  malting  up  the  working  electrode  are  monitored  individually.  These 
techniques  are  often  used  to  probe  the  kinetics  of  mass  transport,  chau'ge  tr<insfer, 
coupled  chemical  reaction  rates,  amd  absorption  (Bard  k  Faulkner,  1980;  R.W.  Mur¬ 
ray,  1986).  In  addition,  AC  impedance  spectra  can  be  recorded  while  the  electrode 
is  pulsed  either  potentiostatically  or  galvanostatically.  These  spectra  are  examined 
in  terms  of  by  an  equivalent  electric  circuit  analog,  thus  indicating  which  process(es) 
control"  the  electrode  behavior  (MacDonald,  1987).  Conway  (1965)  has  represented 
electrode  processes  using  such  analogs.  However,  these  circuits  axe  not  complete  for 
the  hydrogen  absorbing  electrodes  because  the  processes  occurring  in  the  solid  phase 
axe  not  included. 

Also  the  cell  assembly,  shown  in  Figure  4,  can  provide  information  on  the  evo¬ 
lution  of  the  temperature  profile,  T(a;,t),  ajxd  the  concentration  gradients,  Vc(x,t). 
To  measure  T(x,t),  sm«ill  differences  in  temperature  can  be  displayed  by  simply  re¬ 
placing  the  glass  cover  by  a  plastic  sheet  impregnated  with  a  liquid  crystal  (Driscoll 
k  Szpak,  1985).  The  concentration  gradients  in  the  electrolyte  phase  can  be  viewed 
with  the  aid  of  a  Watson  interference  objective.  Using  such  an  apparatus,  one  can 
monitor,  in  real  time,  local  changes  in  the  electrolyte  phase.  An  example,  shown  in 
Figure  5,  demonstrates  the  localized  natxire  of  the  charge  transfer  re2w:tion  and  re¬ 
sultant  concentration  gradients  as  well  as  their  movement  with  the  passage  of  time. 
These  observations,  when  examined  against  the  structured  features  of  the  porous  elec¬ 
trode,  may  lead  to  the  rational  specification  of  the  “effective  parameter,”  i.e.,  in  the 
dependence  between  process  coupling  and  the  geometrical  features  of  the  electrode 
interior. 

4.3.2.  Examination  of  the  Electrode/Electrolyte  Interphase 

The  properties  of  the  electrode/electrolyte  interphase  are  governed  by  solution 
and  electrode  composition  as  well  as  the  metallurgy  of  the  electrode,  e.g.,  grain  size, 
grain  boundaries,  and  the  density  of  defects.  During  charging/discharging,  changes 
in  the  structure  and  appearance  of  the  electrode  surface  occur  and  reflect  the  con¬ 
tributions  of  an  active  electrode/electrolyte  interphaise.  For  example,  Rolison  et  al. 
(Rolison,  O’Grady,  Doyle,  &  Trzaskoma,  1990)  demonstrated  for  Fd  electrodes,  dur¬ 
ing  prolonged  exposure  to  hydrogen  and  deuterium,  that  the  active  participation  of 
the  surface  manifested  itself  as  inhomogenous  ch£inges  in  the  surface  morphology.  A 


Figure  5.  Development  of  concentration  gradients  in  an  experimental  cell,  similu  to 
that  shown  in  Figure  4,  as  displayed  with  the  aid  of  a  Watson  interference  objective. 
Electrode  reaction:  Zn  +2e~;galvemostatic  conditions  j  =  lOmAcm”^; 

distance  X  =  0.3  cm;  the  tim<  between  exposures  is  60  s. 
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realistic  model  must,  therefore,  include  an  active,  changing  electrode/electrolyte  in¬ 
terphase.  This  is  very  difficult  to  accomplish  since  there  are  limited  real-time,  in-situ 
data  on  the  morphological  changes  of  the  electrode  surface  during  charging  and/or 
discharging.  Most  of  the  data  that  is  available  is  from  ex-situ,  post-test  examinatiou 
of  electrodes  which  have  undergone  several  charging/discharging  cycles.  Experiments 
must  be  designed  which  will  provide  the  necessary  data  to  aid  in  developing  the 
model. 

4.3. 2.1.  Relevant  Experimental  Tools.  One  way  to  demonstrate  the  inhomogene¬ 
ity  of  a  surface,  particularly  with  regard  to  absorption,  is  to  view  it  using  Nomarski 
optics  (Szpak,  Mosier-Boss,  Gabriel,  k  Smith,  1994)  so  that  regions  of  preferred 
absorption  can  be  differentiated.  An  example  of  an  experimental  arrangement  is  il¬ 
lustrated  in  Figure  6.  This  arrangement  allows  continuous  in-situ  monitoring  of  the 
electrode  surface  during  charging/discharging.  The  usefulness  of  this  technique  to 
monitor  changes  in  the  Pd/H^O  interphase  during  the  uptake  of  electrochemically 
generated  hydrogen  has  been  demonstrated  (Szpaic,  Mosier-Boss,  Gabriel,  Smith, 
1994). 


Figure  6.  Electrolytic  cell  for  the  examination  of  an  electrode  surface  by  Nomarski 
optics. 

A  -  Pt  counter  electrode;  B  -  Cu  wire  spotwelded  to  Pd  electrode;  C  -  etched  Pd 
surface;  D  -  Mylar  film. 


Since  the  properties  of  the  interphase  are  determined  by  those  of  the  contacting 
phases,  i.e.,  the  electrode  interior  as  well  as  the  electrolyte,  the  extent  to  which  the 
electrode/electrolyte  interphase  facilitates  hydrogen  sorption  can  be  modified  by  chang¬ 
ing  the  composition  of  either  phase.  One  modification  of  interest  is  the  introduc¬ 
tion  of  surface  active  impurities  and  their  effect  on  the  performance  characteristics. 
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Examples  of  additives  which  strongly  modify  the  behavior  of  the  Pdj H^O  system  are 
thiourea,  A^O^.,  and  cyanide  ion.  The  effect  of  these  and  other  additives  on  hydro¬ 
gen  sorption  of  the  AB2  and  AB^  alloys  must  be  ascertained  via  cyclic  voltammetry, 
galvanostatic  and  potentiostatic  pulsing,  and  AC  impedance  techniques.  The  results 
of  these  studies  will  provide  a  greater  understanding  on  how  the  interphase  influences 
hydrogen  sorption  and  will  facilitate  development  of  the  auxiliary  equations  in  the 
model. 

4.3.3.  Experimental  Cell  to  Monitor  Hydrogen  Sorption 

Practical  electrodes  consist  of  particles  of  the  hydrogen-absorbing  material  em¬ 
bedded  in  a  conductive  matrix,  such  as  foamed  Ni.  Data  on  the  hydrogen  absorption 
capacity  and  rates  of  charging/discharging  of  the  material  need  to  be  known  for  com¬ 
plete  modeling.  Such  data  are  also  needed  in  order  to  determine  which  material  is 
best  for  a  given  application  and  to  optimize  the  composition  of  that  material  to  meet 
the  performance  criteria.  These  quantities  for  the  AB2  and  AB^  materials  should  be 
measured  ais  a  function  of  electrode  and  electrolyte  composition. 

The  amount  of  hydrogen  stored  within  the  electrode  structure  and  the  rate  at 
which  it  can  be  transferred  across  the  interphase  in  both  directions  can  be  monitored 
using  a  cell  similar  in  design  as  that  described  by  Will  et  al.  (1991).  This  cell, 
in  modified  form,  is  illustrated  in  Figure  7.  In  this  sealed  cell,  the  MH  working 
electrode  is  separated  from  the  counter  electrode  by  a  fritted  glass  plate  to  prevent 
direct  contact  between  O2  gas  and  H2  gas  generated  during  charging  of  the  catalytic 
walls  of  the  MH  electrode.  After  cell  assembly,  the  space  above  the  electrolyte  is 
filled  will  H2  and  the  pressure  is  recorded.  During  electrolysis  (charging),  H2  gas 
generated  at  the  MH  electrode  either  transports  into  the  solid  or  transfers  into  the 
gas  phase.  The  gaseous  Oj  generated  at  the  positive  electrode  enters  the  gas  phase 
above  the  anolyte  and  is  recombined  on  the  surface  of  a  Pt  catalyst.  The  amount  of 
absorbed  hydrogen  and  the  rate  at  which  absorption  has  occurred  is  calculated  from 
the  pressure  changes  in  the  gas  phase.  To  equalize  the  liquid  level  and  the  pressure 
in  both  compartments,  shunting  is  provided. 


+ 


to  Manometer 


Figure  7.  Design  of  the  experimental  cell.  Tj  and  T2  refer  to  the  temperatures  of  the 
MH  electrode  and  solution,  rt  jpectively. 


5.  CONCLUDING  REMARKS 


The  relevance  of  modeling  to  the  development  of  metal  hydride  technology  can  be 
summarized  as  follows: 

1.  Theoretical  modeling  deals  with  problems  of  general/broad  applicabil¬ 
ity,  i.e.,  it  assumes  a  physical  model  (^rangement)  and  provides  solu¬ 
tions  for  mathematical  statements  (equations).  The  connection  with 
reality  is  thus  limited  by  the  input  parameters,  the  values  of  which  are 
often  unknown. 

2.  To  be  effective  as  a  development  guide,  parameters  relevant  to  the  se¬ 
lected  system  must  be  obtained.  These  parameters,  in  turn,  determine 
the  performance  envelope. 
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3.  In  real  electrode  systems,  mathematical  formulation  may  be  so  complex 
that  unambiguous  interpretation  is  not  possible.  In  such  cases,  exper¬ 
imental  simulation  may  be  helpful  in  displaying  the  most  important 
features,  thereby,  reducing  computation  and  providing  a  more  realistic 
guidance. 

6.  SYMBOLS  AND  CHEMICAL  EQUATIONS 
USED  TO  DERIVE  EXPRESSIONS  IN 

TABLE  I 

6.1.  LIST  OF  SYMBOLS  USED  IN  EQUATIONS  IN  TABLE  I 

a  interfacial  area  per  unit  volume  (cm“^) 

Ca  molar  concentration  of  AlCh  (mole  cm~^) 

c*  molar  concentration  of  electrolyte  (mole  cm“^) 

Ci  molar  concentration  of  species  i  (mole  cm“^) 

c"  initial  electrolyte  concentration  (mole  cm~^) 

Co  molar  concentration  of  solvent  (mole  cm“®) 

Dj  effective  diffusivity  of  SOCh  in  LiCl  film  (cm*  s~^) 

Dia  diffusion  coeflacient  for  acid  (cm*  s"^) 

Die  diffusion  coefficient  for  electrolyte  (cm*  s~^) 

T  Faraday’s  constant  (96,487  C/equiv) 

iok  exchange  current  density  for  reaction  k  (A  cm"*) 

12  superficial  current  density  in  pore  phase  (A  cm"*) 

i„fc  local  current  density  due  to  reaction  k  (A  cm"*) 

I  superficial  current  density  to  an  electrode  (A  cm"*) 
jin  pore  wall  flux  of  species  i  (mole  cm"*  s"') 

ki  rate  constant  of  reaction  (III)  at  anode  (cm  s"*) 

L  electrode  thickness  (cm) 

njt  number  of  electrons  transferred  in  electrode  reaction  k 

Nj  superficial  flux  of  species  i  (mole  cm"*  s~*) 

p  morphology  peirameter 
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R  universal  gas  constant  (8.3142  J  mole"^  K— 1) 

Ro  reaction  rate  of  chemical  reaction  (III)  at  cathode  (mole  cm"^  s— 1) 

s,,o  stoichiometric  coefficient  of  species  i  in  chemical  reaction 

Si^k  stoichiometric  coefficient  of  species  i  in  electrode  reaction  k 

t  time  (s) 

tramsference  number  of  species  i  with  respect  to  the  volume  average 
velocity 

T  absolute  temperature  (K) 

Vi  molar  volume  of  species  i  (mole  cm“^) 

Zi  \  alence  or  charge  number  of  species  i 

Qa  transfer  coefficient  in  anodic  direction 

Oc  transfer  coefficient  in  cathodic  direction 

^LiCi  LiCl  film  thickness  (cm) 

c  porosity  or  void  volume  fraction 

f-LiCi  volume  fraction  of  LiCl 

Tj  ~  4>i  -  <f>2 

r]s  surface  overpotential  (V) 

K  effective  solution  conductivity  (mho  cm~^) 

<r  effective  matrix  conductivity  (mho  cm“^) 

<t>\  electric  potential  in  electrode  matrix  (V) 

<t>2  electric  potential  in  the  solution  (V) 

Subscripts  : 

a  AICI3 

b  LiCl  film/separator  interface 

e  electrolyte 

/  LiCl  film 

k  reaction  k 

p  positive  electrode 

s  separator 
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Superscript  : 


o  initial 

6.2.  REACTIONS  CONSIDERED  IN  DERIVING  EQUATIONS  IN 
TABLE  I 


(I)  Li  Li+  +  c- 

(II)  2Li  +  SOCh  -  2LiCl 

(III)  LiCl  +  A/C/3  Li-^  +  AlCl; 

(IV)  2Li+  +  2c-  +  SOCI2  -  2LiCl  +  ^5  + 

(V)  2A/C/3  +  2c-  +  SOCh  -  2A/C/;  +\S  +  \S02 
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